Friction vs Texture at the Approach of a Granular Avalanche 
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We perform a novel analysis of the granular texture of a granular bed close to stability limit. 
Our analysis is based on a unique criterion of friction mobilisation in a simulated two-dimensional 
packing. In this way, we recover the bimodal character of granular texture, and the coexistence 
of weak and strong phases in the sense of distinct contacts populations. Moreover, we show the 
existence of a well-defined subset of contacts within the weak contact network. These contacts are 
characterized by their important friction, and form a highly coherent population in terms of fabric. 
They play an antagonistic role with respect to force chains. We are thus able to discriminate between 
incoherent contacts and coherent contacts in the weak phase, and to specify the role that the latter 
plays in the destabilisation process. 



Disordered systems exhibit critical behaviours which 
often result in the emergence of rare "catastrophic" 
events for which predictive modelling is problematic. 
Earthquakes, and more generally, rupture of compos- 
ite materials, are well known example of this critical- 
ity 0- Dry cohesionless granular matter, characterized 
by a generic disorder induced by the rigidity of the grains 
and the highly dissipative contacts between them, obeys 
the same type of behaviour. Its response to loading in- 
volves a succession of micro-plastic events which eventu- 
ally lead to a catastrophic failure 0,0]' f° r instance, a 
surface avalanche. 

In spite of their intrinsic disorder, granular materials 
obey a well-defined organisation commonly referred to 
as granular texture Q. In particular, only a fraction 
of the contacts contributes to the quasi-totality of the 
shear stress. This fraction, known as the strong contact 
network and forming force chains, has been extensively 
studied on the ground that it is responsible for the me- 
chanical strength of the whole granular structure. The 
organisation of frictional contacts, i.e. the dissipative 
part of the contact network, has attracted less attention. 
Yet they are responsible for the local dynamics observed 
before failure and are essential for the identification of 
precursors of destabilisation 0, EI • 
In this Letter, we propose an analysis of granular struc- 
ture based solely on the criterion of friction mobilisation. 
Doing so we are able to recover the bi-phasic nature of 
granular texture. More importantly, we show for the 
first time the existence of a well defined population of 
contacts, highly frictional and coherent, and coexisting 
in the weak phase with a large population of incoherent 
contacts. These findings enable us to revisit the charac- 
teristics of the granular texture, and the respective roles 
of the weak and the strong contacts in the destabilisation 
process. 

This analysis was performed by means of two- 
dimensional numerical simulation. We have used a con- 
tact dynamics algorithm [fj assuming perfectly rigid 
grains interacting at contacts through a hard core repul- 
sion and a Coulomb friction law. Beyond the fact that 
contact dynamics treats them as strictly non-smooth, 
these contact laws do not differ from those more com- 



monly used in discrete simulations |l(. 
We consider a granular bed built by random rain of 8000 
circular grains of diameters varying in a ratio 1.5 to in- 
duce geometrical disorder. The solid fraction of the pack- 
ing is 0.82. The height of the bed is 40D and its width is 
200-D, where D is the grain mean diameter. To prevent 
rigid boundary effects, the bed is periodic in the horizon- 
tal direction, namely no vertical walls are introduced to 
confine the packing. A constant rotation rate is applied 
(corresponding to a rotation of 10 -3 deg per time step), 
gradually bringing the slope 9 of the bed from initially 
zero to the limit value 9 C for which a surface avalanche is 
triggered. 

The Coulomb friction law imposes an upper limit on the 
tangential force fx at contact. This upper limit (the 
Coulomb threshold) is given by /i//v , where /at is the nor- 
mal force at contact and \x is the coefficient of friction. A 
necessary condition for slip motion to occur between two 
grains in contact is fx = M/-/V- ln the- following, the value 
of /i is set to 0.5 and its influence, though not investigated 
here, will be discussed later. A coefficient of restitution 
allows for the modelling of energy exchanges during bi- 
nary collisions; since we are interested in micro-plasticity 
only, we set its value to zero and consider purely dissipa- 
tive contact interactions. 

The mobilisation of friction at each contact is simply 
measured by the ratio of the tangential to the normal 
force times the coefficient of friction: 



with n G [0, 1] following the friction law. The case rj = 1 
coincides with the Coulomb threshold and a probable 
slip motion at the contact. We compute the mean value 
(77) over the total number of contacts N c in the packing 
as a function of the slope angle 9 (Fig. QJi). The mean 
normalised kinetic energy of the grains as a function of 
9 is displayed in Fig. [Tfr. The onset of the avalanche 
occurs for 9 = 9 C ~ 20°. While the mean level of 
friction increases regularly, we observe some dynamical 
activity intensifying few degrees before 9 C . This inter- 
val corresponds to a metastable state during which lo- 
cal dynamical rearrangements occur more frequently due 
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FIG. 1: a) Mean mobilisation of friction forces (??)jv c aver- 
aged over all the contacts and b) mean kinetic energy Ek 
normalised by the typical potential energy of a grain of mass 
m for a height D as a function of the slope angle 9. 

to an enhanced mobilisation of friction 0,0. In the 
following, we thus focus on this interval of slope angles 
= [6 C — 3° , 9 C ] , and the quantities discussed hereafter 
are averaged over it. 

We consider now different populations of contacts fol- 
lowing their friction level r\. We introduce a parameter £ 
to which the friction level r\ a is compared for each contact 
a. A contact a belongs to if i] a > (1 — £): 



C s = {a GN c /r, a > (!-£)}. 



(2) 



The parameter £ varies in the interval [0, 1]. The number 
of contacts belonging to is N$. The proportion p^ = 
N^/N c of contacts belonging to is plotted against £ 
in Fig. |2 For £ close to (i.e. for a high mobilisation of 
friction), p^ reaches 18%, exceeding by far the case of a 
uniform distribution. This value bespeaks the existence 
of a local peak in the distribution of r\. Contacts such 
that £ — > thus form a well denned subset. In the same 
way, the increase of p^ when £ — > 1 is more rapid than 
a uniform distribution would allow, showing again the 
non-uniformity of the distribution of friction at contacts. 
We have checked that these features were robust with 
respect to /i and the boundary conditions in piles close 
to stability limit. 

A picture of the mobilisation of friction depending on 
contact orientation is obtained when plotting the value 
of r\ averaged over angular sectors if + Atp for contacts 
directions: 
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(3) 



a=r 



where N£ is the number of contacts for which the ori- 
entation of the unit vector n normal to the contact plan 
belongs to [ip,<p + Ay>]. The function r] v (tp) is plotted 
in Fig. |3 in polar coordinates. The slope of the granular 
bed coincides with the horizontal direction, and the direc- 
tion of gravity is represented by a dashed line (at —70°). 
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FIG. 2: Proportion of contacts belonging to the subset of 
contacts Cj. The dashed line represents the case of a uniform 
distribution of r\. 




FIG. 3: Polar representation of the mean friction mobilisation 
r) v as a function of the orientation of the contacts. 



We observe a higher mobilisation of friction for contacts 
directed at ~ 25° [%], i.e. nearly perpendicular to the 
gravity. Note that this angular map of friction mobilisa- 
tion does not need to compare with the angular map of 
contact orientations. 

To compare the anisotropy of friction mobilisation with 
the anisotropy of the geometrical contact network, we 
compute partial fabric tensors defined over the con- 
tact subsets Cf p~0l | : 



$i = — 



i 



a=l 



(4) 



where n a is the unit vector of contact a, and i and j refer 
to coordinates x and y (with x parallel to the free sur- 
face, and y perpendicular to the free surface and pointing 
downward). The major principal direction 9^, of the ten- 
sor corresponds to the mean direction of the contacts 
in regarding the slope. The difference of the eigen- 
values gives the anisotropy of the contact sub-network 
Q. These two quantities are displayed in Fig. 0] against 
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FIG. 4: a) Anisotropy of the contact network in the sub- 
set and b) major fabric direction 0$ (counted positively 
anticlockwise) as a function of £ (see text). The dashed line 
shows the major stress direction over all the contacts. 

£. We see that the largest anisotropy occurs for £ ~ 0, 
i.e. for high friction mobilisation. It then declines reg- 
ularly and passes by its lowest value (close to zero) for 
£ ~ 0.7 before increasing again up to £ ~ 1. This transi- 
tion coincides with a rotation of 7r/2 of the fabric tensor. 
We can thus discriminate between three groups of con- 
tacts on the unique criterion of friction mobilisation: 

1. Contacts with a very high mobilisation of fric- 
tion (£ — y 0), i.e. contacts where slip motion is 
likely to happen. Their orientations obey a maxi- 
mum anisotropy (0.75), three times superior to the 
anisotropy of the whole packing, and define a mean 
direction at ~ 40°, orthogonal to the major stress 
diection (~ —50°). These contacts are thus forming 
a highly coherent subset in terms of fabric. 

2. Contact with intermediate level of friction mobil- 
isation (0.1 < £ < 0.7): they very weakly con- 
tribute to the geometrical anisotropy. Accordingly 
their mean orientation allows only for a small devi- 
ation from the direction dictated by the first group 
of contacts. Their small contribution to the fabric 
suggests that they form an incoherent ensemble. 

3. Contact with a very low mobilisation of friction 
(£ — * 1): they strongly contribute to the geometri- 
cal anisotropy. Accordingly, their direction allows 
for a dramatic rotation of the fabric tensor, forc- 
ing the latter to follow the major stress direction. 
Their strong contribution to the fabric show that 
they form a coherent ensemble. They represent 40% 
of the total number of contacts (see Fig. [2J. 

Discriminating between contacts on the ground of 
contact friction discloses features strongly reminiscent 
of the bimodal character of granular texture p|. The 
granular texture is characterized when discriminating 
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FIG. 5: a) Mean force transmitted by contacts of C% as 
a function of £ and b) Probability density function Po.05 of 
F/Fh within the subset of contacts Co.os- 



between contacts on the ground of the intensity of the 
force they transmit. In the case of our granular bed, 
the effect of gravity has to be filtered out to distinguish 
between weak and strong contacts. This is done by 
normalizing all contact forces F by the mean contact 
force Fh at the depth h corresponding to the position 
of the contact. The mean value F^ of the normalised 
force transmitted by the contacts of is plotted in 
Fig. as a function of £ : higher mobilisation of friction 
force (small £) coincides with weaker forces and vice 
versa. This concordance between force transmission 
and friction mobilisation implies that the transitions 
observed in Fig. 0] partly reflect the transition from the 
weak to the strong contact network. 
Focusing on contacts for which the mobilisation of 
friction is the highest, we consider the subset C0.05, is. 
the contacts satisfying n > 95%. For these contacts, 
we plot in Fig. |SJd the distribution of the intensity 
Po. 05 (.F/Fh) °f the forces they transmit. We observe 
that most contacts of Co.os transmit very weak forces, 
such that F/{F) h < 0.25. In any case, 90% of the 
contacts of Co.os belong to the weak contact network, 
namely transmit forces such that F/Fh < 1< 

The interplay between force transmission and friction 
mobilisation is illustrated in Fig. where one snapshot 
of the state of the granular pile at an instant of the 
interval Ag — [8 C — 3°,# c ] is shown. In gray are rep- 
resented the force chains defined by the strong contact 
network, with line width varying proportionally to the 
force intensity. Black segments of fixed width represent 
the contacts belonging to Co.os- The following features 
are apparent: contacts of Co.os ar e belonging to the 
weak phase (namely none of them is part of a force 
chain), they show a strong geometrical anisotropy (for 
we distinguish indeed a leading orientation), and their 
direction is orthogonal to the main stress orientation (i.e. 
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FIG. 6: Sample of the granular bed in the interval [9 C — 3°, 9 C ]- 
In gray are represented the force chains (contacts such that 
F/Fh > 1) and in black contacts such that n > 0.95/i. 



contacts of C0.05 are in their great majority transverse 
to force chains). 

On the unique criterion of friction mobilisation (as de- 
fined by equation lfl]l) in a packing close to avalanching, 
we are able to recover the bimodal character of granu- 
lar packing and the existence of two distinct phases, in 
agreement with the analysis in terms of a strong and a 
weak contact network. Moreover, we show the existence 
of a well-defined population of contacts within the weak 
contact network. These contacts are defined by their very 
high mobilisation of friction. They form a highly coher- 
ent population, showing a strong anisotropy and a well 
defined mean direction. In this sense, they strongly dif- 
fer from other weak contacts. Their orientation, orthog- 
onal to force chains, and their susceptibility to slip due 
to their level of friction, impart to them an antagonist 
role with respect to force chains. Frictional contacts and 
strong contacts form the extremal state of our analysis, 
and they play extremal roles in terms of stability. By con- 



trast, all other contacts in the weak network are highly 
incoherent. 

This analysis relies on the relative value of the normal 
and tangential force at contacts, and is independent of 
any mean quantity or any cut-off value reflecting the 
mean state of the packing. As such, it purely reflects 
the intrinsic organisation of the packing. However, the 
level of friction mobilisation is necessarily measured with 
respect to the coefficient of friction /1. In particular, the 
proportion of highly frictional contacts is expected to de- 
crease when the value of /1 increases. An investigation of 
the effect of this parameter on the results reported here 
should be undertaken. However, in the same way that 
the characteristics of force transmission are robust with 
respect to the material properties, we expect the organ- 
isation of friction to be a robust feature of granular tex- 
ture. 

Back to the description in terms of strong and weak con- 
tact networks, the analysis of granular texture in terms 
of friction allows to revisit the role of these two entities. 
The contribution of the strong contact network to the 
stress state designates it as solid skeleton of the granu- 
lar structure, and responsible for its stability. By con- 
trast, the weak contact network appears as an "inter- 
stitial" phase essentially screened from external loading, 
and proping the force chains. In this work, we show that 
the weak contact network is actually hosting a subset 
of contacts whose role in destabilisation is fundamental. 
Although forming a less coherent phase when analysed 
on the ground of force transmission (i.e. representing 
a small anisotropy), the weak contact network contains 
an intrinsic coherence when analysed on the ground of 
friction. This makes the weak contact network an active 
agent of the destabilisation process. In other words, the 
macroscopic behaviour of the granular packing is partly 
dependent on a phenomenology which escapes the mean 
description in terms of stress state. 
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Fellowship FP6 program. 



[1] M. P. Allen and D. J. Tildesley, Computer Simulation of 
Liquids, Oxford University Press, 1997 

[2] J. V. Andersen, D. Sornette and K. Leung, Phys. Rev. 
Lett. 78 11, 2140 (1997) 

[3] A. Daerr and S. Douady, Nature 399, pp. 241-146 (1999) 

[4] F. Darve and F. Laouafa, Numerical Models in Geoma- 
terials, Balkema, pp. 85-90, 1999 

[5] S. Deboeuf, E. M. Bertin, E. Lajeunesse and O. Dauchot, 
Eur. Phys. J. B 36, p. 105, 2003. 

[6] M. Jean and J. -J. Moreau, in Proceedings of Contact Me- 
chanics International Symposium PPUR, 1992 

[7] A. Kabla, G. Debregeas, J.-M. Di Meglio and T. Senden, 



cond-mat/0412074 
[8] F. Radjai, D. E. Wolf, M. Jean, and J.-J. Moreau, Phys. 

Rev. Lett. 90, 61 (1998). 
[9] F. Radjai, H. Troadec and S. Roux, Granular Materials: 
Fundamentals and Applications, RS.C, 2004. 

[10] M. Satake, IUTAM symposium on deformation and fail- 
ure of granular materials, Balkema, 1982. 

[11] L. Staron, J. -P. Vilotte, and F. Radjai, Phys. Rev. Lett. 
89, 204302 (2002). 

[12] D. Volfson, L. S. Tsimring, and I. S. Aranson Phys. Rev. 
E 69, 031302 (2004) 



